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Abstract 

The total cross section of the top quark-antiquark pair production near threshold 
in 77 collision is computed analytically up to the next-to-leading order in perturbative 
and nonrelativistic expansion for general photon helicity. The approximation includes 
the first order corrections in the strong coupling constant and the heavy quark velocity 
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1 Introduction. 



Theoretical study of the top quark-antiquark pair production near the two-particle threshold 
is based on the key observation that the relatively large width T t of the top quark serves as a 
natural infrared cutoff for long distance strong interaction effects 0. Quantitatively, the top 
quark width is mainly saturated by the decay channel t — ► Wb and is equal to T t = 1.43 GeV. 
Near the production threshold the nonrelativistic approximation is accurate and the relevant 
scale y/r t mt = 15.8 GeV, m t being the top quark mass mt = 175 GeV, is much larger 
than Aqcd that makes perturbative QCD applicable for the theoretical description of the 
threshold top quark physics 0, On the other hand the Next Linear Collider provides an 
opportunity of experimental study of high energy 77 interactions which can be used for the 
top quark production. Using the laser backscattering method one can obtain 77 colliding 
beams with the energy and luminosity comparable to those in e + e~ collisions |4j]. Because 
of high precision which can be achieved in such experiments and the high accuracy of the 
theoretical description, the process of top quark pair production in 77 collisions near the 
two-particle threshold has been recognized as a promising place for thorough investigation of 
quark interactions and for determining the numerical values of the strong coupling constant 
a s , the top quark mass m t , and the top quark width T t . Moreover, the strong dependence 
of the cross section on the photon helicities extends possibilities of studying the details of 
the top quark threshold dynamics [0]. 

As is well known the ordinary finite order perturbation theory of QCD breaks down in the 
threshold region of particle production and the resummation of the Coulomb effects is neces- 
sary that can be systematically done in the framework of nonrelativistic QCD (NRQCD) ||. 
The threshold effects are treated exactly by using the nonrelativistic Schrodinger equation 
with effective parameters. This equation can be also obtained from the corresponding Bethe- 
Salpeter equation in the nonrelativistic limit. Relativistic effects for the top quark are rather 
small and can be treated as corrections. Note that the characteristic scale of the Coulomb 
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effects for the top production a s m t is comparable numerically with the cutoff scale y/T t m t 
so these effects are not suppressed by the top quark width. Quantitatively, Coulomb effects 
require to sum up power terms in the variable (up to color group factors) a s \Jm t /T t ~ 1.1 
which is not small and the coresponding function is not accurately approximated by the 
finite Taylor series. The determination of the higher order strong coupling corrections and 
relativistic corrections in this case requires the perturbative expansion to be performed near 
the Coulomb approximation rather than near the Green functions of a free theory. Techni- 
cally, it makes the calculation rather involved. Only recently an essential progress has been 
achieved in analytical evaluation of the next-to-leading order (NLO) and the next-to-next- 
to-leading order (NNLO) corrections to the heavy quark vacuum polarization function near 
the threshold which enters the analysis of the heavy quark threshold production in e + e~ 
annihilation || [7], || |], |nj [H| . The corrections are found to be relatively large in the case of 
the top quark production and very large in the case of the bottom quark threshold sum rules. 
One can expect the similar effect in the 77 threshold production of heavy quarks. Therefore, 
from phenomenological point of view, the calculation of a s and relativistic corrections to the 
corresponding cross section is necessary for accurate quantitative study of the process. A 
part of the NLO corrections (the one to the cross section for the colliding photons of the 
same helicity) has been studied in ref. H using the numerical algorithm [Q]. The complete 
analytical result for the corrections is still absent. 

In the present paper we formulate the relevant (NRQCD) framework for the calculation 
of the higher order corrections and obtain analytically the NLO corrections which include 
the first order corrections in the strong coupling constant and the heavy quark velocity to 
the nonrelativistic Coulomb approximation for general photon helicity. The analytical result 
in the case of the opposite helicity colliding photons is of special importance because the 
available numerical methods are not applicable in this case ||. We evaluate the relativistic 
T t a s corrections to the cross section for the opposite helicity colliding photons that cannot be 
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found within the pure nonrelativistic approximation and require an additional information 
on the relativistic structure of the full theory. The procedure of matching to the full theory is 
used to recover the necessary parameters of effective NRQCD. We present also the analytical 
NNLO result for the parameters of the resonances in the cross sections for 77 — > tt and 
e + e~ — > tt top quark threshold production. 

The paper is organized as follows. In the next Section we formulate the nonrelativistic 
approximation for the total cross section of top quark pair production in 77 collisions. In 
Section 3 corrections to the nonrelativistic Green function and its derivative which enter 
the expression for the cross section are considered. In Section 4 we discuss the effects of 
nonvanishing top quark width. In Section 5 we present the NNLO analysis of the resonance 
structure of the cross section. The last Section contains the results of numerical analysis and 
our conclusions. 

2 The nonrelativistic expansion for the cross section. 

We fix our notation of kinematics for the top quark pair production in 77 collisions for 
working out the nonrelativistic expansion for the cross section as follows. The collision 
process is given by 

7(9i > e(9i)hfe, e(q 2 )) -> f (pi, <7i)i(p 2 , <r 2 ) 

where qi,q 2 are momenta of initial photons, pi,p 2 are momenta of the final state fermions 
(top quark-antiquark pair), qi + q 2 = pi + p 2 , (<?i + ffc) 2 = s is the square of the total energy 
of two photon beams in the photons zero-momentum frame. Here e A (gj) are the photon 
polarization vectors, <7i )2 are spin variables of the fermions. After reducing the photon states 
the amplitude of 77 — > tt transition reads 

M(77 -> tt) = {eQ t )h»{ qi y{q 2 )(tt\T^) (1) 
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where eQ t , Qt = 2/3, is the top quark electric charge. The quantity in eq. ([I]) is a Fourier 
transform of the T-product of two top quark electromagnetic (vector) currents J\ = t^\t 

T, u (q) = i J TJ^x/2)J u {-x/2)e^dx , (2) 

where q = (qi — q2)/2 is the difference between the phonon momenta. 

The T-product of the electromagnetic currents in eq. (§) can be expanded into a series 
over the local operators using a consistent expansion in x at small x within the Wilson 
operator product expansion |~3|]. Being inserted in eq. (0) it converts into a power series in 
the variable l/{m 2 — q 2 ). Because q 2 = — s/4 where s is a total photon energy squared the 
expansion is not singular near the fermion-antifermion threshold where s ~ 4m 2 . In fact, 
in the threshold region q 2 ~ — m\ and the corresponding expansion becomes an expansion 
in l/mf. The leading term of the operator product expansion (a contribution of the local 



operator of dimension three) reads |13 



T ( 3) = 2ie^ a/3 g a t7 /3 7 5 t 



^ m 2 -q 2 ' (3) 



Next contributions to the transition amplitude are related to the dimension four local oper- 
ators of the operator product expansion which have the following explicit form 



T (4a) = Zt 7(M Du) t . ■ 

V m 2 -q 2 ' 



rp{4b) _ Da t . s 

^ ~ (m 2 t -q 2 ) 2 ' ' 



r ( 4c ) = i2q a tq (jl y u) D a t 



^ (m 2 -q 2 ) 2 ^ 

where D a is a covariant derivative containing the gluon field. Note that operator (^]) does 
not contribute to the amplitude because q^e^qi) = 0. The operator product expansion is 
organized in such a way that near the two-particle threshold the matrix element of a higher 
dimension operator between the vacuum and ti state is suppressed by either the quark pair 
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energy E = \fs — 2m t counted from the threshold or explicitly by powers of a s . Thus, for the 
transition amplitude the operator product expansion converts into a simultaneous expansion 



in a s and in (3 = y 1 — 4mf / s, the parameter of nonrelativistic near-threshold expansion. The 
operators in eqs. (0-|6|), however, are relativistic ones. Therefore, to obtain the consistent 
nonrelativistic expansion of the amplitude one has to expand these operators in a series 
over the nonrelativistic operators in the framework of NRQCD [14]]. The nonrelativistic 



expansion of the operators entering eqs. (|| |4j) which are necessary for the further analysis 
up to a perturbative normalization factor reads 

t 7 V = <P*X + ■ ■ ■ (7) 
itj(i Dj) t = i^a {i Dj)X + ■■■ (8) 

where <p (x) is the nonrelativistic two component top (antitop) quark spinor, o"j is the Pauli 
matrix and ellipsis stands for the higher dimension operators. Note that only the time 
component of the current in eq. (^[) and space components of the tensor in eq. (|J) contribute 
in the leading order of nonrelativistic expansion. 

Let us turn to the calculation of the cross section. We deal with the normalized cross 
section 

R(s) a(77 ^ 



a(e + e~ — ► ) 
where the standard lepton cross section 

/ + - + -n Ana QED 

a{e^e -» /i> ) = — 

3s 

is just a kinematic normalization factor. 

In the case of the same helicity colliding photons only the operator of dimension three 
contributes to the amplitude up to NLO. Using the leading term of the nonrelativistic ex- 
pansion for the quark current in eq. (0) one arrives at the following representation for the 
cross section [H 



R ++ (E) = 6QtN c C+ + (a s ) ( ^ ) ImG(0, 0, k) (9) 



where G(x, y, k) (k 2 = —m t E) is the nonrelativistic Green function which accounts for the 
soft (threshold) effects and 



4 J n 

is the NLO perturbative coefficient matching correlators of relativistic and nonrelativistic 
currents jnjj with a s being taken at "hard" normalization scale \Xh- The higher dimension 
operators in the expansion of eq. (|2]) and in the nonrelativistic expansion (|7|) lead to relative 
0(f3 2 ) corrections to the cross section. 

For the opposite helicity colliding photons the amplitude up to NLO is determined by 
the operator given in eq. @. By using the leading term of eq. (§) for the corresponding 
cross section one gets [@ 

a x a y ImG(x,y,A;)|, = o,, = o (10) 

with the perturbative coefficient being [ITB 



C h (a s ) = l-A . 

7T 

Note that the leading contribution to R + ~(E) comes from the quark pair with the orbital 
momentum I = 1 because I = state cannot be produced by the opposite helicity photons 
with the total angular momentum J = 2 while the dimension three operator (|3|) interpolates 
only I = states and, therefore, does not contribute to R + ~ cross section. The dimension 
four operator (|5|) does not contribute to R + ~ in the leading order because it has the parity 
P = — 1 while the quark pair with the orbital momentum I = 1 in the final state has the 
parity P — 1. In full analogy with the R ++ case the higher dimension operators lead only 
to relative 0{[3 2 ) corrections to R + ~ cross section. 

The nonrelativistic Green function appearing in eqs. (|9|JT0|) in the considered approxima- 
tion satisfies the Schrodinger equation 

(-— + V c {x) + A.Vix) + —) G(x, y, k) = *(x - y) (11) 



where x = |x|, Vc(x) = — Cpa s (fi s )/x is the Coulomb potential with fi s being the soft 
normalization scale which determines the strong coupling constant. In the NLO one has 
to take into account the first order perturbative QCD corrections to the Coulomb potential 
Ai V which is of the following form |L9| 

A 1 V(x) = ^Vc(x)(Cl + C\ H^s)) = ^V c (x)Cl ln(x Ml ) (12) 

where 

Cl =a x + 2(3 o7E , C\ = 20o, 
ai = y Ca ~ y Tirn /' ^ = ^ ex P ifll C i) ■ 

The color symmetry SU(3) group invariants for QCD are Ca = 3, Cp = 4/3, Tp = 1/2, 
nf = 5 is the number of light fermion flavors, and /3q = 11Ca/3 — 4Tpnf/3 is the first 
coefficient of /3-function. Here ^ E = 0.577216 ... is the Euler constant. 

Note that up to an overall normalization factor the expression for R ++ up to NLO 
coincides with the normalized cross section of the top quark threshold production in e + e~ 
annihilation. 

In zero order in a s eqs. (PJiCf) are reduced to Born approximation for corresponding cross 
sections and contain only continuous spectrum starting from the two-particle threshold 

R ++ (f3) = 6QtN c /39(s - Am 2 t ) + 0((3 3 ) } 

R+~(p) = 8QtN c f3 3 6{s - Amf) + 0{(3 5 ). 

The full relativistic result for the cross sections R ++ {(3) and R + ~{j3) in tree approximation 
which sums up all powers of (3 was obtained in ref . [[I7j] . At the one loop level the analytical 
result has been recently obtained as a power expansion in (3 up to (3 10 order |18| . 
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3 The nonrelativistic Green function beyond the lead- 
ing order. 

In this section we describe the technique for computing the nonrelativistic Green function 
beyond the leading order. The corrections are found within the perturbation theory around 
the Coulomb solution as a leading order approximation relevant for computation of the cross 
section near the threshold. The developed technique is applicable for computing correc- 
tions to the nonrelativistic Green function with any orbital momentum I and is a direct 
generalization of the method suggested in \$, j9j pTXf] . 

The nonrelativistic Green function has the standard partial wave decomposition 

oo 

G(x,y,k) = £(2l + l)(xy) l P l (xy/xy)G l (x 1 y 1 k) (13) 
1=0 

where Pi(z) is a Legendre polynomial. The partial waves of the Green function of the pure 
Coulomb Schrodinger equation G c (x.,y,k) are 

nC( r „ k\ - mtk (ok\ 2l P - k ^+y) V L% l +l (2kx)L% L +l (2ky)m\ 

where v = X/k, X = a s Cpm t /2, L^(z) is a Laguerre polynomial 

e z z~ a ( d \ m 
ml \dz J 

and a s is taken at the soft scale fi s . At y = the partial waves take the form 

Gf(x, 0, k) = —{2k) 2l e- kx T{l + 1 - u)U(l + 1 - v, 21 + 2, 2kx) (15) 
2ix 

where U(a, b, z) is the confluent hypergeometric function. 

Only I = component of the Green function contributes to its value at the origin 

G(0,0,fc) = G o (0,0,Jfe). 

For the Coulomb Green function we find an explicit expression in the form 

G c (x, 0, fc)Uo = ^ - 2Aln (2xn f ) " 2A ^ + ^ 



+2 7B - 1 + *i (1 - u))) (16) 

where ^f n (z) = d n In Y(z)/dz n and T{z) is the Euler T-function. In the short distance limit 
x — > the Coulomb Green function G c (x, 0, k) has 1/x and lnx divergent terms. These 
terms, however, are energy independent and do not contribute to the cross section. Hence 
they can be subtracted. In eq. (|T6| ) fif is an auxiliary parameter. Though an accurate 
definition of this parameter is important in the NNLO analysis J7|, [H| it does not enter the 



cross section of top production for the stable top quark, or in zero width approximation. 
The derivative of the Green function at the origin is saturated with its I — 1 component 

d x <9 y G(x, y, k)\ y= o tX =o = 9Gx(0, 0, k) 

and for the Coulomb Green function we have 

+A ( 2( p _ A2) + ln ( A) + 27E _ £ + (1 _ „,) + ^ . (17) 



One can immediately see that this expression has no lowest pole present in eq. (p!6|). In the 
short distance limit x — > the derivative of Coulomb Green function (or partial wave with 
/ = 1) has l/x n (n = 1,2,3) and lnx singularities. In contrast to Coulomb Green function 
itself some of the divergent terms now are k dependent but also do not contribute to the cross 
section for T t = because they have no discontinuity across the physical cut in the complex 
energy plane in the zero-width approximation. The case of non-zero top quark width needs 
more detailed analysis given in the next section. 

The solution to eq. ([□]) with the logarithmic correction to the Coulomb potential can be 
found within the standard nonrelativistic perturbation theory around the Coulomb Green 
function taken leading order approximation 

G(x, y, k) = G c (x, y, k) + A 1 C(x, y, k), 
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A 1 G(x, y,k) = -J G c (x, z, k) A 1 V(z)G c (z, y, k)dz . (18) 

The calculation goes along the lines developed in refs. || [| |TTJ5 Only I = component of 
eq. (|HJ) is necessary for the calculation of the correction to the Green function at the origin. 
The integral in eq. (|IBD diverges at x, y = 0. This divergence is induced by the singular 1/x 
term in eq. (|T6|) which is related to the singular behavior of the free (a s = 0) Green function 
G F (x, y, k) in the limit y — 0, x — > 0. The divergent part can be separated in the following 
way. Eq. (|18D for x, y = can be written in the form 

A 1 G(0, 0, k) = A 1 G (0, 0, k) = - y (G%(0, x, k) - Gq (0, x, k)) 2 A 1 V{x)d^ 

+2 J (G%(0, x, k) - G^ (0, x, k))G% (x, 0, k)A l V(x)dx + J G% (0, x, k) 2 A 1 V{x)d^j . (19) 

By using representation (|H|) for Gq(x, y, k) and the same expression with a s = for 
Gq(x, y, k) one gets 

2 



AiG (0, 0, k) = - ^ £ F(m)F(n) J e~ 2kx L 1 m (2kx)L 1 n (2kx)A l V(x)d^ 

«|/W /f^MMM.x, /^f^*) (20) 



m=0 

where 

F(m 



(m + 1) (m + 1 — u) 
and we used the properties of the Laguerre polynomial 

r(m + q + l) ^ L« ( Z ) 

L - (0) = r( a + i)r (m + i)' J? ^ = * r(a) - 

Only the last term is divergent in eq. (p0|) . After regularization the divergent integral in that 
term becomes 

roo e -2kx } n (^ x ) X 



Jo x 2 



1 Th.e correction to the Green function at the origin has been also studied using the numerical algorithm 
in refs. § | g |. 
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where L(k) = — ln(2/c//ii) and ellipsis stands for the divergent part. The divergent part 
being k independent does not contribute to the spectral density. Two finite integrals in 
eq. (|19D are 

roo (m + l)^i(m + 2), m = n 

/ e~ ZL n(z) L ln( Z ) H Z )zdz = I n + 1 

J0 , m> n 

V m — n 

roo 

/ e- z L] n (z)\n{z)dz = _2 7s -* 1 (m+l). (21) 
Jo 

To compute the first integral we rewrite it in the form 

jf e- z L\{z)Ll{z) Hz)zdz = ±(J e' z L^L^z^dz) ^ . (22) 
By using the relations 

L'(z)-T r ^- tt + n ) L * ( Z ) 
^ T{/3 - a)T{n + l) 

f°° -zTaf \Tcct \ a, x T(m + a + l) 

/ e z L a n {z)L a m (z)z a dz = 5 mn K 

Jo 1 [m + 1) 

for (3 = 1, a = l + e the integration in the right hand side of eq. ( p2[) can be performed 
analytically. Then taking the derivative in e at e = we get the first line of eq. (|2~ip . The 
second integral in eq. ( |21| ) can be computed using the same technique. Thus the final result 
for the correction is 

AxGoCO, 0, fc) = E *W V + 1) (L(k) + 2)) - 2 E E F H 

Z7F ^ \m=0 m=l n=0 

XjP ( n )^^_ + 2 £ F(m) - 2 7i; - *i(m + 1)) - Tjj L(A;) + -L{k) 2 . (23) 

Note that the sum in the finite limits in this equation can be reexpressed through some com- 
bination of ip functions of different arguments. This however is not important for numerical 
analysis at present when the actual calculation is performed with computers. 
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The calculation of the correction to the derivative of the Green function at the origin (or 
partial wave with I = 1) is a bit more cumbersome. The corresponding expression reads 

a x a y A 1 G(x, y, k)\ Xiy=Q = 9AiGi (0,0, k) = -9 J Gf(0,x,k) 2 A 1 V(x)x 2 dx. (24) 

Eq. (|T7|) contains l/x n (n = 1, 2) singular terms with the coefficients depending on a s which 
lead to the divergence of the integral in eq. (p4j). Hence this is not enough to subtract just 



the free Green function from the Coulomb one to separate the divergent part of eq. ( 24]) as 
has been done in the analysis of the correction to the Green function at the origin. The 
analog of eq. dl9|) now reads 



where 



AxG^O, 0, k) = - (y (Gf {x, 0, k) - D(x, k)) 2 A l V(x)x 2 d^ 
-2 J (Gf (x, 0, k) - D(x, k))D(x, k)A 1 V(x)x 2 dx + J D(x, k) 2 A 1 V(x)x 2 dxj (25) 



D(x,k) = ^^ e ~ kx H{2kx), 



By using the relation 



m=0 

eq. fl2"5"|) can be transformed into 



2(m+l) m + 2 2(m + 3) 



AiGi(0,0,fc) = - ( — — ) V F(m)F(n) ( e- 2kx L 3 n (2kx)L 3 n (2kx)A 1 V(x)x 2 dx 



+2 £ ^(™) / e~ 2te L^(2£;x)#(2£;x) A^x^x + / e - 2kx H(2kx) 2 A 1 V(x)x 2 dx j (27) 



m=0 

where 

z^ 2 - i: 



F(m) 



(m + 2 - u){m + l)(m + 2)(m + 3) 
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and the divergence is contained in the last term. After a regularization this term can be 
directly computed with the result 



I(k) = - ^ Je-^H(2k X fA lV ( X ) X ^ = - (7 *~ 1)2 - g 



- (4 - 3 lE )v + i^±M±l!, 2 + Iz±V + i^V (28) 

+ (7* - 1 - 31/ + ^i^V + ^ 3 + i, 4 ) + (-- + ^ 2 ) L(fc) 2 + . . . 

where dots stand for divergent term which has no imaginary part and does not contribute 
to the spectral density in zero width approximation. A set of finite integrals appearing in 
eq. fl25|) can be computed by the technique which has been already applied for the analysis 
of the correction to the I = partial wave. The results for relevant integrals read 

poo I (m + + 2)(m + 3)^i(m + 4), m = n 

J Q e- z Ll(z)L 3 m (z)\n(z)z 3 dz= I (n + l)(w + 2)(w + 3) 

\ m — n 

POO 

J (m) = / e- z L 3 m (z) \n(z)dz = -2*i(m + 1) - 4 7e + 3, 

JO 

roo 

Ji(m)= e- z L 3 m (z)\n(z)zdz = (m + l)(-*i(m + 2) - 2 7J5 + 2), (29) 
Jo 

Mm) = / o °° e-*L*Jz) \n(z)z 2 dz = {m + ^ + 2) (-^(m + 3) - 2 7i; + |) . 

Substituting these formulae in eq. (p^ ) we arrive at the final result for the correction to the 
I = 1 partial wave at the origin 

A 1 G 1 (0, 0, k) = a f° X ™f { V F(m) 2 (m + l)(m + 2)(m + 3) (L(fe) + *x(m + 4)) 

Z7T 187T 1 



\m=0 



-2 £ I] 1 F(m)F(n) ( n + 1 )( n + 2 )( n + 3 ) + 2 g ^ (m) ( 2 J (m) + (m + l)(m+ 2)L(k) 



m=l n=0 m 72 m=0 



+ (1 + i/)(Ji(m) + (m + l)L(k)) + V[V ^ l \ j 2 (m) + 2L(fc))) + . (30) 

Note that it is quite straightforward to extend our method to calculation of the correction 
to a component of the Green function with arbitrary I. 
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The Green function in every order of the nonrelativistic perturbation theory has to be 
written in the form which includes only single poles in the energy variable and is more natural 
for the Green function of a nonrelativistic Schrodinger equation. For the Green function at 
the origin this representation reads 

Gm ) ~Lw^ + -Jo -w^w dE (31) 

where ipm^'ify is the wave function at the origin, E$ is a I = bound state energy, the sum 
goes over bound states and the integration in the second term is performed over the state 
of continuous part of the spectrum. In this way the corrections to Green function stemming 
from the discrete part of the spectrum can be reduced to the correction to the Coulomb 
bound state energy levels 

^ = -^TTTF + S 2C >' [i(m) + + 2)1 ) (32) 



and to the values of Coulomb bound state wave functions at the origin 

'«°>' 2 = ^TW i 1 + £ CJ ( 3i < m > " 1 " 2 ^ + ^TT + *'(™ + 2 > 

-2(m+l)* 2 (m + l))) (33) 

where 

L(m) = In ( . 

\ C F a s m b J 

For the derivative of the Green function at the origin one has 

d x d y A lG (x, Y, E) \ x , y =o - L gl) _ E _ iQ + -J E >-E-iO dE (34) 

where 

\^'m,E'(0)\ 2 = d^ !E ,(x.)d y 1p m! E>(y)\x,y=0 

and E$ is a I — 1 bound state energy. In NLO approximation these quantities read 

E - - -^{&W (' + £ 2c; [Um + + * l(m + 4)l ) • (35) 
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and 



A 5 (m + l)(m + 3) / a 



|C(0)| 2 = + [ l + [~\ {m + 2 ) - 1 + 5L (™ + !) + 5 ^(™ + 4 ) 

^ (n + l)(n + 2)(n + 3) \\ 

~ (m + l)(m + 3)(m-n) 2 jj ' 1 j 

The continuum contributions in eqs. ( |3~T| , ^) can be directly found by subtracting the discrete 

part of these equations expanded around the Coulomb approximation up to NLO from the 

result obtained within the nonrelativistic perturbation theory for the Green function and its 

derivative at the origin. 
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4 Effects of the finite width of top quark. 

As has been already mentioned only sufficiently large t-quark decay width suppresses the 
nonperturbative effects of strong interactions at large (~ 1/ Aqcd) distances and makes the 
perturbation theory applicable for the description of the t-quark threshold dynamics. The 
accurate description is achieved in the framework of NRQCD that provides the adequate 
leading order approximation for the problem. However main features of the physical situation 
can be reflected in a simpler approach to which we restrict ourselves in the present paper. 
According to the prescription of ref. |l| the non-zero width can be taken into account by 
direct replacing E — > E + iT t in the Green function near the production threshold. This is 
justified because the threshold dynamics is nonrelativistic one and is rather insensitive to 
the hard momentum details of t-quark decays. This prescription is not equivalent to the full 
NRQCD description in higher orders of perturbation theory in the strong coupling constant. 
Also, strictly speaking this prescription is valid in its simple form only when the complete 
factorization of the hard and soft contributions to the cross section takes place. This is 
realized in the case of R ++ cross section. In the case of _R H cross section, however, the 
nonrelativistic approximation is not able to describe properly the entire effect of the non- 
zero top quark width [0] . Indeed, eq. (|H]) for a nonvanishing width in the limit x — > has the 
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divergent imaginary part with the leading power singularity ~ Ft fx related to the free Green 
function singularity and logarithmic singularity ~ T t a s \nx produced by the one Coulomb 
photon exchange. This clearly indicates that the coefficient in front of the constant term 
linear in T t in the cross section gets a contribution from large momentum region and cannot 
be obtained within the nonrelativistic theory. This constant can be obtained within the 
full relativistic theory. The correct treatment of the problem consists in matching NRQCD 
expression for this part of the cross section to the one obtained in full theory by direct 
inserting the (complex) mass operator into the t-quark propagator. In the leading order in 
a s this has been done semiphenomenologically in refs. [fj, 0. In the present paper we limit 
ourselves to formulas obtained in 0, [L2| that are sufficient from phenomenological point of 
view. The contribution linear in the top quark width reads 

> 2 

mi 

= 8QfN c T t p a ((9m 2 -p 2 )(m 2 -p 2 )f 2 T t (p 2 )dp 2 = Q ^ 8QfN c T t 
m t J(m w +m b ) 2 {rn? t - p 2 ) 2 + 4m 2 T t (p 2 ) 2 32nT t mf ' m t 

where p 2 is the Born amplitude squared, dp is the relativistic phase space of two unstable 

particles for E = (see ref. |T^]), T t (p 2 ) is the imaginary part of the t-quark mass operator 

(r t (m 2 ) = T t ) and the leading t — > Wh decay channel is taken into account. In eq. ( p7|) 

m\y = 80.3 GeV is the M^-boson mass |2(J and mj = 4.8 GeV is the 6-quark pole mass [|TTf . 

In the final expression for the cross section this contribution has to be multiplied by the 

hard normalization coefficient C^~(a s ). Numerically this contribution to the cross section 

is suppressed in comparison with one of the regular pure nonrelativistic terms of eq. ( |i~7D 

which saturate the total result. 

To find the coefficient in front of the T t a s term in the imaginary part of eq. (|i~7D or, 

equivalently, to fix the parameter /i/ one has to know the one-gluon correction to Born 



2 

R+- = 64irQfN c ( ^dp 
' J mt 



approximation. This correction can be obtained from QED result [21 1 and for E = is equal 
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to 7rA|p|. A corresponding contribution to the cross section reads 

A|p| 



J 1 1 L+ 



8QtN c \T t ( 9m t ~ V 2 ){< - P 2 ) r t (P 2 W _ 8QfN c \T t /m o Q _ 



raj J(m w +m b ) 2 (mf — p 2 ) 2 + 4m 2 r t (p 2 ) 2 8r 4 m 2 m 2 V V T 



(38) 

where the term with the logarithmic dependence on T t is separated and computed analytically 
within the approximation T t (p 2 ) = T t . Following the general line of the effective field theory 
approach one has to equate the contribution to the cross section determined by the 0(T t a s ) 
terms in the imaginary part of eq. flT7| ) at E — to the relativistic expression (|3"8|). This 
equation results in the following matching relation 



which fixes the auxiliary parameter /i/ = 0.13m t . Note that in this relation the logarithmic 
dependence of eq. (|T7| ) on Y t exactly matches one of eq. (|38| ) so /i/ is related only to the hard 
scale of the process as one expects from the general ground. 

The above analysis cannot be considered as a completely relativistic one - Born amplitude 
and the Coulomb one-gluon correction are taken in the leading order of the nonrelativistic 
expansion though the relativistic phase volume is used. This, however, is justified because 



the integrals in eqs. (|37], |38|) are saturated within the region |p| ~ y/m t T t «m t (p ~ mf) 
where the nonrelativistic approximation works well. By the same reason the function T t (p 2 ) 
in these equations can be well approximated by its constant value T t . 

Note that for non- vanishing width some of the singular terms in eq. (p8|) which describes 
the correction to the Green function also have imaginary part and contribute to the cross 
section. This contribution, however, has very weak energy dependence and is suppressed by 
an extra power of a s . Thus, this is beyond the accuracy of our approximation. 
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5 Resonances of the cross section in the next-to-next- 
to-leading order. 

Though the complete NNLO analysis of the cross section is not still available, some important 
parameters of the cross section can be found in NNLO approximation. Indeed, the position of 
the poles of the nonrelativistic Green function which correspond to resonances in the cross 
section can be determined within the nonrelativistic approximation P, [TO, O, This 



has been explicitly done for the spin triplet states in refs. ||T0| , |TT| , P2| . This analysis can be 
directly generalized to the spin singlet states which are produced in the two photon collisions. 
The only difference is in the correction induced by the part of the Breit potential responsible 
for the hyperfme splitting. We, however, do not need the complete result because the most of 
resonances are smoothed out by the relatively large top quark width. The numerical analysis 
shows that only the ground state resonance in I = partial wave cross section R ++ (E) is 
distinguishable. Its separation from others is not completely covered by the infrared cutoff 
provided by the top quark width. Indeed, using the pure Coulomb formulas for estimates 
within the order of magnitude we find 

\E^-E[ 0) \ = 1^0.6 GeV 

to be compared with the top quark width T t = 1.43 GeV. The spacing between next radial 
excitations for I = partial wave (which is also equal to the first one for I = 1 partial wave) 
is much smaller 

\E[ 0) - 4 0) | = - E?\ = -— « 0.11 GeV 

36 m t 

and is completely smeared out with the top quark width. Therefore it is numerically justified 
to treat separately only the first resonance in I — partial wave cross section and to sum 
all corrections to the denominator of the Green function to avoid the appearance of the 
double pole in the correction which is large at the resonance energy. For all other states this 
procedure makes no numerical difference because the infrared cutoff is sufficiently large to 
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smear out all picks in the cross section. For R + ~(E) cross section the lowest contributing 
state has I = 1 and is rather close to other states and therefore is indistinguishable as a 
separate contribution. 

Using the results of refs. || [10, 11], it is straightforward to find the energy of this 
resonance in NNLO approximation counted from the threshold 

E res (ll - tt) = = -— (l + ^2C\ (L(A) + 1 — Je) + ^ (2C 2 1 (L 1 (\) + 1 - lE ) 

m t V Air 47r V 

Heir - lE ) 2 + 1 - y - * s (i)) + 2C\ ((l 2 (\) + 1 - lE f - 1 + y) ) 

+CWs(^ + ^]). (39) 



C F 16/ 



where 



Co = ( ^ + 47l J 01 + 2(A + 2/3oai) 7£ + a a , 



v 3 

C\ = 2(A + 2A,a a ) + Cl = 4/3 2 , 

/4343 „ 2 tt 4 22 w „,\^o /1798 56 



02 = \ w2 + 6?r x + t c(3) r^Ur^ y c(3) J ^ 

y-iecis^c^T^+^T^) 2 , 

and £(z) is the Riemann ^-function. This value is related to the energy of the bound state 
resonance of the top quark production in e + e~ annihilation by the hyperfine splitting 

4 A 2 

E res {e + e- -> tt) - E res (ll -> tt) = tj — C|a 2 . (40) 

This NNLO result provides us with the information on the convergence of the pertur- 
bative expansion and, in particular, restricts the allowed region of the normalization point 
which can be used for reliable estimates. In fact, the normalization points of a s entering 
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the coefficients Ch and the nonrelativistic Green function can be different when NNLO cor- 
rections are considered. The difference between the normalization points of the hard and 
soft corrections can be noticed only in higher orders of perturbative expansion. This gives 
an additional possibility to improve the convergence of the perturbation theory. The typical 
hard scale of the problem is the heavy quark mass m t . The soft physical scale of the prob- 
lem is determined by the natural infrared cutoff related to the top quark width y/m t T t that 
measures the distance to the nearest singularity in the complex energy plane and/or by the 
characteristic scale of the Coulomb problem A i.e. fi s ~ 15 GeV. For top quark both scales 
are rather close to each other that makes possible a uniform description both perturbative 
QCD and Coulomb resonance effects. Indeed, for fi s ~ 15 GeV the NLO correction to the 
resonance energy reaches its minimal magnitude (see Fig. 1). However, from Fig. 1 we see 
that at this scale the NNLO correction is large and the series for the resonance energy di- 
verges. At first glance this seems to contradict our physical intuition. However, since the 
normalization scale is defined in a rather artificial MS scheme that has little to do with pecu- 
liarities of ti physics being originally designed for describing a massless quark approximation 
in e + e~ annihilation, there is no reason for a literal coincidence of fi s parameter with any 
physical scale of the process. The relative weight of the NNLO corrections is stabilized at 
fi s ~ 40 GeV where the series for the resonance energy looks like E res = E^(l + 0.6 + 0.5). 
Moreover, for the larger values of the normalization point the /i dependence of the cross 
section becomes small (see the numerical analysis of the next Section). Thus our analysis 
implies the optimal choice of normalization point in the range /i s ,>40 GeV. The convergence 
of the series, however, is not fast and the uncertainty related to the truncation of the se- 
ries exceeds those due to the nonperturbative effects which are suppressed parametrically 
as (Aqc£)/A) 4 p3| . This implies relatively large NNLO corrections to the R ++ cross section 
(similar effect takes place in e + e~ — > ti top quark threshold production 0, §])• Note that 
the apparent convergence for the resonance energy can be slightly improved by changing the 
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definition of the top quark mass and using MS scheme for defining the mass parameter of 
the top quark instead of its (perturbative) pole mass the convergence of the series for the 
resonance energy becomes better. 



6 Summary and conclusion. 

Thus in the present paper the total cross section of the top quark pair production near 
the threshold in 77 collision is computed analytically up to the next-to-leading order in 
perturbative and nonrelativistic expansion for general photon helicity with the top quark 
width being taken into account. 

To demonstrate the significance of the NLO corrections we plot the functions R ++ (E) 
and R + ~(E) in Born (a s = 0), leading orderQ and NLO approximations in Fig. 2 and Fig. 3. 
We use m t = 175 GeV, r t = 1.43 GeV and a s {mz) = 0.118 as typical numerical values of 



corresponding parameters pOfl , fih = mt for the hard normalization scale and /i s = 25 GeV, 
50 GeV, and 75 GeV for the soft normalization scale (in the leading order approximation 
the upper curves correspond to ji s = 25 GeV and lower curves correspond to /i s = 75 GeV). 
The results are almost independent of [ih and their dependence on the soft normalization 
scale decreases with increasing 

The main results that we have obtained from the study of the NLO corrections to the 
total cross section of the top quark-antiquark pair production near threshold in 77 collision 
are: 

• The typical size of the NLO corrections is 20% for R ++ (E) and 10% for R + ~(E). 

• Inclusion of the NLO corrections leads to considerable stabilization of the theoretical 
results for the cross sections against changing the normalization point. 



2 In the leading order approximation + = = 1 and the Coulomb Green function is used. 
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• Relatively large top quark width smears out all the resonances in R + ~(E) and only 
the ground state resonance survives in R ++ (E). The inclusion of the higher order 
corrections makes this peak more distinguishable. Thorough experimental study of 
the location and shape of this pick can provide precise data for extraction of the mass 
and width of top quark and the numerical value of the strong coupling constant. 

• The NNLO correction to the resonance energy in 77 — > tt top quark threshold produc- 
tion is relatively large that implies large NNLO corrections to the R ++ cross section. 
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Figure captions 

Fig. 1. The relative weight of the NLO corrections E^°/E^° - 1 (solid line) and NNLO 
corrections E^ LO /E^° - 1 (dotted line) to the resonance energy as a function of the 
normalization point /x s . 

Fig. 2. The normalized total cross section R ++ (E) of top quark pair production for the 
colliding photons of the same helicity in leading order (dotted lines) and NLO (bold solid 
lines) for /i s = 25 GeV, 50 GeV and 75 GeV. The normal solid line corresponds to the Born 
approximation. 

Fig. 3. The normalized total cross section R + ~(E) of top quark pair production for 
the colliding photons of the opposite helicities in leading order (dotted lines) for /i s = 
25 GeV, 50 GeV and 75 GeV and NLO (bold solid lines) for fi 8 = 50 GeV (the NLO 
result is given for only one normalization point because it has very weak dependence on the 
normalization scale). The normal solid line corresponds to the Born approximation. 
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